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Abstract 

It was demonstrated recently in Bychkov et al., Phys. Rev. Lett. 101 (2008) 164501, that the physical mechanism of flame 
acceleration in channels with obstacles is qualitatively different from the classical Shelkin mechanism. The new mechanism is 
much stronger, and is independent of the Reynolds number. The present study provides details of the theory and numerical 
modeling of the flame acceleration. It is shown theoretically and computationally that flame acceleration progresses noticeably 
faster in the axisymmetric cylindrical geometry as compared to the planar one, and that the acceleration rate reduces with 
increasing initial Mach number and thereby the gas compressibility. Furthermore, the velocity of the accelerating flame saturates 
to a constant value that is supersonic with respect to the wall. The saturation state can be correlated to the Chapman- Jouguet 
deflagration as well as the fast flames observed in experiments. The possibility of transition from deflagration to detonation in 
the obstructed channels is demonstrated. 



1. INTRODUCTION 

In the process of deflagration-to-detonation transition 
(DDT) in tubes with a closed end, a slow premixed flame 
accelerates spontaneously from the closed end and trig- 
gers detonation [THT3]. A qualitative explanation of the 
process was first proposed by Shelkin, Ref. [1 j. Accord- 
ing to this scenario, the burned gas expands and drives 
a flow in the fuel mixture. The flow becomes nonuni- 
form because of the no-slip boundary condition at the 
wall which, together with turbulence, distorts the flame 
front, increases the burning rate, and leads to the accel- 
eration. An accelerating flame front pushes compression 
waves that continuously heat the fuel mixture ahead of 
it until an explosion is triggered that eventually develops 
into a detonation. 

While the Shelkin mechanism relies prominently on the 
action of turbulence, it was recently shown theoretically 
that flame acceleration is possible even in its absence, 
in tubes with smooth adiabatic wall p~4j [15] . This the- 
ory has been validated by extensive numerical simula- 
tions in Refs. [I4j [15] and supported by experiments 
in smooth micro-tubes [16]. The theory and modeling 
further demonstrate that laminar flame acceleration be- 
comes quite weak in wide tubes with increasing Reynolds 
number of the flow, and as such loss to the wall may ac- 
tually terminate the process. Thus obstacles placed in 
the tube [TJ I5HTU1 [T2] appear to be an essential factor in 
order to overcome the loss and consequently support the 
DDT. It is generally believed that these obstacles gener- 
ate stronger turbulence, which increases the burning rate 
and facilitates the flame acceleration. 

However, in our recent work [17J we have demonstrated 
that the obstacles play a more important role than just 
producing turbulence, in that they provide a specific 
physical mechanism of flame acceleration that is quali- 



tatively different from the Shelkin mechanism. This new 
mechanism is extra strong, providing flame acceleration 
that is independent of the Reynolds number, through the 
tube width, and as such may be quite important for tech- 
nical applications. Specifically, flame propagation in an 
obstructed channel creates pockets of fresh fuel mixture 
between the obstacles, as shown in Fig. [T] Gas expansion 
due to delayed burning in the pockets produces a power- 
ful jet flow in the unobstructed part of the channel. The 
jet flow renders the flame tip to propagate much faster, 
which produces new pockets, generates a positive feed- 
back between the flame and the flow, and leads to flame 
acceleration. The accelerating flame reaches supersonic 
speed with respect to the tube wall, and consequently 
triggers explosion and detonation. 



The present paper extends the work of Bychkov et al., 
Ref. [17 , in which the basic concepts of the new mech- 
anism were outlined, by presenting details of the theory 
and simulations. Specifically, we discuss the influence 
of the planar /axisymmetric flow geometry and gas com- 
pression on the acceleration rate. We show that the flame 
accelerates much faster in cylindrical tubes than in planar 
channels, and that the acceleration is slowed down due to 
gas compression. Furthermore, as the Mach number of 
the flow increases, the acceleration process saturates to 
statistically steady flame propagation at supersonic speed 
with respect to the tube wall. This saturation state may 
develop prior to the attainment of explosion and then 
detonation. The flame speed in this state may be corre- 
lated with the Chapman- Jouguet (CJ) deflagration speed 
[T8H20] and with the state of fast flames observed exper- 
imentally [6j [2T] , Finally, we demonstrate numerically 
that flame acceleration may lead to explosion and deto- 
nation triggering. 




FIG. 1: Flame propagation in a channel with obstacles (a) 
and the jet generation mechanism (b). 



2. THEORY OF FLAME ACCELERATION 

Figure [I] is a schematic of the problem under study. 
The flame propagates from the closed end of a semi- 
infinite channel of half-width (radius) i?, with a frac- 
tion a < 1 blocked by obstacles. The central part of 
the channel, of half- width (1 — a)R, is unobstructed. 
The flame propagates extremely fast along the unob- 
structed part of the channel, leaving behind pockets of 
unburned mixture, between the obstacles, which will be 
burned later. The deep narrow spaces between the obsta- 
cles, namely the pockets, act as mini-channels in which 
the flame can be considered to propagate mainly in the 
radial direction. This assumption is most appropriate 
when the obstacles are placed close to each other with 
deep pockets, Az « aR, and with slip boundary con- 
ditions at the wall. In general, burning in the pockets de- 
pends on the obstacle geometry, whose influence will be 
discussed later. We nevertheless recognize qualitatively 
the same mechanism of flame acceleration in simulations 
and experiments involving complicated obstacle shapes 

GanniiE]. 

We now briefly state the fundamental concepts of the 
new mechanism presented in Ref. [17 . We employ 
the standard model of an infinitesimally thin flame front 
propagating normally to itself with the unstretched lam- 
inar flame velocity Uf. The tube wall is ideally slip and 
adiabatic. Non-slip boundary condition is not needed 
for flame acceleration in the new mechanism and the 
Reynolds number is not involved in the calculations. At 
the initial stage of flame acceleration, the flow may be 
treated as incompressible, 



V • u = 0, 



(1) 



while the fresh gas trapped in the pockets is burning. 
Expansion of the burnt gas is characterized by the density 
ratio of the fuel mixture and the burnt gas, = pf/pb, 
which is quite large for most flames, with = 5 — 8. In 



the model of a planar laminar flame front in the pockets, 
the fuel mixture in a pocket is at rest, while the burnt 
gas is pushed out with the velocity (9 — l)Uf, see Fig. [I] 
(b). This value determines the gas velocity at the border 
of the unobstructed part of the channel 



(G-l)Uf at 



±(l-a)R. 



(2) 



The shape of the flame tip is of minor importance for the 
present mechanism, and it may be taken to be planar at 
all times. Consequently, the flow of the burnt gas in the 
unobstructed channel part is potential. Accounting for 
boundary condition we find the flow velocity 



K; u z ) = — r^r {-x] z). 

(1 — a)R 



(3) 



The propagation speed of the infinitesimally thin flame 
front with respect to the burnt gas is 



dZ 



f 



dt 



u z (z f ) = eu f . 



Using the velocity distribution (J3|, we find 
dZ f _ (9 - l)U f 



dt 



(l-a)R 



■z f + eu f . 



(4) 



(5) 



Integrating Eq. ^ with the initial condition Zf(0) = 0, 
we obtain a strong exponential acceleration of the flame 
tip 

(T^ = e^i[ex P W/i?)-i], (6) 

with the scaled acceleration rate 



9 



1 



a 



(7) 



The derivation of Eqs. Q - ^ explains the basis of the 
new acceleration mechanism presented in Ref. [17] . This 
mechanism is quite powerful, with the acceleration rate 
([7]) much larger than that of the Shelkin mechanism in 
smooth tubes [14 . The scaled acceleration rate does not 
depend on the Reynolds number, and hence the viscosity 
and the tube width. As will be discussed below, viscos- 
ity and turbulence may induce supplementary effects in 
the new acceleration mechanism. As pointed out in Ref. 
[T7] , this new acceleration mechanism has many features 
in common with the acceleration of finger flames, Refs. 
[22j [23] , although the finger flame acceleration is quite 
limited in time, yielding maximum increase in the burn- 
ing rate by a factor of only 10-15 relative to the planar 
flame speed [23]. In contrast, the new physical mecha- 
nism leads to supersonic flame propagation with respect 
to the tube wall, with possible transition to explosion 
and detonation. It was also demonstrated in Ref. [T7] 
that the new mechanism remains effective even when the 
flame skirt touches the main wall of the channel. 
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FIG. 2: Typical obstacle shapes discussed in ICDERS-2007. 



Theoretical understanding of the new acceleration 
mechanism enables the analysis of the optimal obstacle 
geometry for flame acceleration. Figure [2] shows typical 
obstacle shapes discussed in Refs. pUT3]. Specifically, in 
Fig. [2] (a) we have rectangular obstacles of thickness d 
and spacing Az. In this case delayed burning occurs only 
between the obstacles, while the volume occupied by an 
obstacle itself is ineffective for flame acceleration. Aver- 
aging volume production in the burning process over the 
obstacle step, we rewrite the boundary condition Q as 

\(u x )\ = (l-d/Az)(e-l)U f at x = ±(l-a)R, (8) 



which leads to the reduced acceleration rate 
a = (l-d/Az) Q ~ 1 



I- a 



(9) 



Thus, thinner obstacles produce stronger acceleration. 
Another possible factor suggested was obstacle phase 
shift at the bottom of the channel in comparison to that 
at the top. This suggestion, however, is not expected 
to substantially affect the laminar flame dynamics, al- 
though it could influence the turbulence generated in the 
flow in later stages. The triangular obstacle shape in Fig. 
[2] (b) also renders the flame acceleration slower because 
it reduces the volume of the fresh fuel mixture trapped 
between the obstacles. Consequently, the optimal design 
for flame acceleration is that of the infinitely thin obsta- 
cles shown in Fig. [I] For a fixed tube radius R, maxi- 
mum acceleration rate is achieved for maximum possible 
blockage ratio, i.e. for minimum value of (1 — a)R. How- 
ever, for a fixed width of the free channel part (1 — a)R, 
the acceleration rate does not depend on the total tube 
radius R, i.e. it does not depend on the depth of the 
pockets. This is true, of course, only when pockets are 
sufficiently deep with a comparable to unity. Very small 
obstacles (a <C 1) cannot be treated as obstacles and 
the new mechanism is ineffective. Still, typical experi- 
mental configurations employ obstacles with the block- 
age ratio of 1/4 < a < 3/4, for which the present mecha- 
nism is quite effective. The dimensional acceleration rate 
(0 — l)Uf / R(l — a) is determined by the expansion factor 
0, the laminar flame velocity Uf and the half- width of 
the unobstructed channel part R(l — a). The theoretical 



model of a laminar flow does not predict any dependence 
of the acceleration rate on the spacing between the obsta- 
cles. As we shall demonstrate numerically in Sec. 4, the 
spacing between the obstacles determines primarily the 
amplitude of the velocity pulsations (see also Ref. [T7]). 
Still, these pulsations do not change the average acceler- 
ation rate, which remains quite close to the predictions 
of the theoretical model. 

Additional increase in the acceleration rate is realized 
in the axisymmetric geometry. Here we consider an ax- 
isymmetric tube with obstacles in the form of planar 
rings blocking the space (1 — a)R < r < R similar to 
Fig. [I] Suppose the flame tip accelerates as Zf — Zf(t). 
A pocket between the obstacles at the position z starts 
burning at the instant tf(z), where tf(z) is the inverted 
function Zf(t). The flame in the axisymmetric pockets 
expands with the radius growing as 



R f = (l-a)R+Uf[t-t f (z)]. 



(10) 



The radial velocity at the exit of the pocket, at r = (1 
a)R, is given for an incompressible flow as 



(0 - l)R f U f = -(1 - a)Ru r , 



(11) 



which determines the boundary condition at the border 
of the unobstructed part of the channel r = (1 — a)R, 

Ur = _ (e - i) Uf (1 + j^^[t - */(*)]) . (12) 

In the event of exponential, or near-exponential, flame 
acceleration 



Z f = Z [exp(aU f t/R) - 1], 



(13) 



the instant tf(z) varies logarithmically, and hence slowly, 
with z, as 



R 



]n(*/Z + l), 



(14) 



where Zq is some amplitude. 

Let us first determine the flame acceleration neglecting 



the increase in the flame radius Rf of Eq. (10) in com- 
parison with (1 — a)R. Then, the axisymmetric solution 
to the continuity equation ([!]) in the unobstructed part 
of the channel is 



(15) 



which is the axisymmetric counterpart of Eq. ([3|. The 
solution Q - ([7]) is respectively modified in the axisym- 
metric geometry as 



dZ f (Q-l)U f 

^T =2 (l-a)R Zf+GUf > 



(16) 



3 



with 



z f = e 

(l-a)R 2(6-1) 



and 



[exp(aU f t/R) - 1], 



9-1 



1 



(17) 



(18) 



Equation (18) shows that flame acceleration in the ax- 
isymmetric geometry is twice that in the planar case. 

We next account for the increase in i?/ as compared 
to (1 — a)R. Taking the radial velocity component in the 
same form as in Eq. ([l5]), 



(e - i)u f 

u r = -^rz -r^-r 



1 



U 



f 



(l-a)R V (l-a)R 
we find the respective z-velocity component 

^ _je-i)u f 



[t - t f (z)] , (19) 



z + 



(l-a)R V ( 
and its value at the flame tip 



dz 



(l-a)R 



(l-a)R 



<*-*/(*)> 



where 



(t - t f (z)} = t 



\- J tf( z ) dz - 



(20) 
(21) 

(22) 



When the fuel mixture in the first pocket is almost 



averaging in Eq. (22) yileds 

Of. 



completely burnt, accounting for Eqs. (14), (17), the 

Zo 



1 



■In 



Z 



1 



(23) 



where <jq is given by Eq. (18). The second term in Eq. 



(23) diminishes asymptotically to zero with time, but the 



first term leads to corrections to Eqs. (16), (18) as 

1 



dz f = (e-i)u f 

dt (1 - a)R 



2(6-1) 



and 



a = 2 



6-1 



1 



a 



2(9-1) 



Z f + eU f , (24) 



(25) 



The second term in Eq. ( |25| ) may be treated as correc- 
tions to the first one in the limit of 2(6 — 1) >> 1, which 
holds with good accuracy of about 7% for realistic fuel 
mixtures. Furthermore, this correction becomes impor- 
tant only when the leading part of the flame skirt has 
almost touched the main wall, when burning in the first 
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FIG. 3: Resolution test: flame tip velocity versus tip po- 
sition for Ma = 0.001, a = 1/2 and various mesh sizes 
(0.125; 0.25; 0.5)L/. 



pocket is almost finished. At the beginning of the burn- 
ing process, the acceleration rate should be approximated 



by Eq. (18). However, even during this initial period of 
flame acceleration, the flame velocity may attain fairly 
large values relative to the sound speed, at which the 
incompressibility assumption fails. Thus, the flame ac- 



celeration rate should be evaluated by Eq. (18) rather 
than Eq. ([25]). 

The theory above constitutes the backbone of the 
new mechanism of extremely fast flame acceleration in 
tubes/channels with obstacles. Still, this theory may be 
developed further to incorporate other effects, such as gas 
compression, viscosity, non-slip at the wall, etc, which 
also influence the acceleration process. Here, we discuss 
briefly some of these effects. 

(1) Gas compression: The present theory is developed 
for incompressible flows. Compressibility is expected to 
moderate the flame acceleration considerably by recog- 
nizing that the flame velocity is limited by the CJ defla- 
gration value. Preliminary numerical results on this sub- 
ject are given in Ref. [17], and similar effects for smooth 
tubes are also discussed in Ref. [20]. We shall perform 
additional modeling later to demonstrate the role of gas 
compression. 

(2) Curved flame shape in the pockets and at non-slip 
wall: Flame shape in the pockets is not necessarily pla- 
nar, as we assumed in the calculations. Inclined or curved 
flame shapes provide faster propagation of the front and 
stronger cumulative expansion. A flame may acquire a 
curved shape because of intrinsic instabilities, viscous ef- 
fects, etc. For example, the viscosity and non-slip wall 
can increase the flame speed by a factor of about 1.5 [24J. 
The curved flame shape in the pockets renders accelera- 
tion of the flame tip stronger. In such cases the planar 
flame velocity in Eq. ^ should be multiplied by some 
factor describing the increase in the flame velocity in the 
pockets. 

(3) Turbulence in the pockets: Flow in the unob- 
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FIG. 4: Snapshots of temperature and velocity in the flow 
generated by an accelerating flame in the planar geometry for 
the following parameters: (a, b) Ma = 0.001, Az/R = 1/4, 
a = 2/3 (appears in [T7]); (c, d) Ma = 0.005, Az/R = 2, 
a = 1/3; (e, f) Ma = 0.005, Az/R = 2, a = 2/3; (g, h) 
Ma = 0.001, Az/R = 1/4, a = 1/2 (appears in [17]). 



structed part of the channel pushed by the accelerating 
flame also generates turbulence in the pockets, see, for 
example, [6j [10l [T7J. Turbulence corrugates the flame 
front and increases the burning rate in the pockets. The 
effect of turbulence becomes especially noticeable when 
the distance between the obstacles is comparable to or 
larger than the obstacle size, Az ex aR. 

(4) Large spacing between the obstacles: Large spacing 
Az causes two competing effects. On the one hand, the 
flow in the unobstructed part of the channel becomes less 
confined, which may be interpreted as decreased effective 
blockage ratio and thereby reduces the acceleration rate. 
On the other hand, large spacing leads to stronger turbu- 
lence in the pockets thus augmenting flame acceleration. 



Numerical simulation [17 showed that these effects com- 
pensate each other at moderate values of Az/R =1, 2, 
with turbulence having a slightly stronger effect in the 
developed stage of flame acceleration. 

(5) Loss to the wall: In smooth tubes heat loss to the 
wall reduces expansion of the burning gas dramatically 
and slows down flame acceleration. However, in tubes 
with sufficiently deep and thin obstacles this negative 
influence may be reduced considerably. Indeed, thin ob- 
stacles may be heated rapidly and do not extract much 
energy from the flow. At the same time, loss to the main 
wall is minimized since it is mostly in contact with the 
fresh fuel mixture instead of the burnt gas. One should 
expect that loss to the main wall is equivalent to an ef- 
fective smaller radius of the tube, R e ff < R, but with 
the same size of the unobstructed part of the channel 
(1 — a)R. Since the acceleration rate is determined by 
(1 — a)R, not by i?, then one should expect minimal in- 
fluence of the loss in such a configuration. 



3. BASIC EQUATIONS AND NUMERICAL 
METHOD 

We performed direct numerical simulation of the hy- 
drodynamic and combustion equations including trans- 
port processes and Arrhenius kinetics. Both 2D pla- 
nar and axisymmetric cylindrical flows were investigated 
using the Navier- Stokes system of the governing equa- 
tions presented, e.g., in Refs. [HI [15]. To avoid the 
thermal- diffusive instability we took unity Lewis num- 
ber Le = 1, with Pr = 0.75 and the dynamical viscosity 
ji = 1.7 x 10 _5 Ns/m 2 . The fuel-air mixture and burnt 
gas were assumed to be a perfect gas with a constant 
molar mass m = 2.9 x 10 _2 kg/mol. We considered a 
one-step irreversible Arrhenius reaction with an activa- 
tion energy E a , pre-exponential factor of inverse time 
dimension r^ 1 , first order dependence on concentration 
of the fuel mixture, and first or second order dependence 
on density similar to Ref. [26 . The first order was used 
in all studies of flame acceleration and saturation to the 
CJ state. The second order was employed to obtain ex- 
plosion and transition to detonation, since in that case 
flame dynamics is much more sensitive to pressure and 
temperature build-up in the compression waves gener- 
ated by an accelerating flame. In the simulations we took 
E/R p Tf = 32 in order to have better resolution of the 
reaction zone. In the present study we focus mainly on 
the flame acceleration and preheating of the fuel mix- 
ture ahead of the flame front, which do not depend on 
the activation energy. The activation energy is crucial for 
time and position of explosion triggering and DDT. How- 
ever, proper quantitative investigation of DDT cannot 
be performed within a simplified one-step mechanism of 
chemical reaction kinetics. Instead, it requires separate 
description of the chemical kinetics at high temperatures 
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FIG. 5: Profile of the scaled z- velocity along the channel axis 
plotted for a = 2/3, Ma = 0.001, Az/R = 1/4 at the time 
instants Uft/R — 0.02 — 0.16 equally spaced in time. The 
theoretical line is related to Eq. Q. 



(flame, detonation) and low temperatures (explosion ig- 
nition). For this reason, in the present work we perform 
only qualitative study of explosion triggering choosing a 
scaled activation energy convenient for such a study. The 
factor tr was adjusted to obtain a particular value of the 
planar flame velocity Uf by solving the associated eigen- 
value problem [271 [28]. For example, taking the planar 
flame velocity Uf = 34.7 cm/s, we set tr = 4.06 • 10 -8 s 
for the first-order reaction. The flame thickness is con- 
ventionally defined as 



L 



f 



PrpfUf' 



(26) 



where pf = 1.16 kg/m 3 is the unburnt mixture den- 
sity. We took initial temperature of the fuel mixture 
Tf = 300if , initial pressure Pf = 10 5 Pa, adiabatic in- 
dex 7 = 1.4, and initial expansion ratio = 8. We 
took different values of the initial Mach number in the 
range 10 -3 < Ma = Uf/c s < 10 -2 , with the lower limit 
corresponding to realistic methane and propane flames. 
By varying the Mach number, we investigated the in- 
fluence of gas compression on flame acceleration. The 
theory of Sec. 2 does not involve the Reynolds number, 
which implies minor dependence of the results on the 
channel/tube width, provided it is sufficiently large. On 
the other hand, simulation of burning in wide channels 
with obstacles can be quite time consuming. For this rea- 
son, in the simulations we used a moderate value of the 
channel half- width (radius), R = 2ALf. Some of the 2D 
planar simulation runs were performed only for half of 
the channel assuming symmetry. However, the majority 
of runs were for the entire channel, which is especially 
important for the developed stages of flame acceleration 
involving turbulence as well as for explosion triggering. 

The present modeling is also relevant to combustion 
in micro-channels, which is a rapidly developing subject. 
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FIG. 6: Profile of the scaled x- velocity along the cross-section 
z/R = 0.5 plotted for a = 2/3, Ma = 0.001, Az/R = 1/4 at 
the time instant Uft/R — 0.16. The theoretical line is related 
to Eq. (TO. 



Recent theory, modeling and experiments [T4HT6] demon- 
strated the possibility of DDT in smooth micro-tubes. 
The present results also successfully predict flame accel- 
eration and DDT in micro- channels with obstacles, while 
recognizing that the present analysis of course has a wider 
implementation than micro-channel flows. According to 
the theory of Sec. 2, the same extra- strong mechanism of 
flame acceleration works even in very wide tubes, though 
it is impossible to attain such a state in numerical simula- 
tions at present because of the inevitable computational 
limitations. The channel width employed in the present 
simulation determines the Reynolds number related to 
the laminar flame speed Re = UfR/v = R/LfPr = 32. 
The Reynolds number related to the flow Re = (u z )2R/u 
can be larger by several orders of magnitude due to flame 
acceleration and thermal expansion of the burnt gas. 
Flame acceleration and increase in the Reynolds num- 
ber may produce turbulence in the flow. Indeed, our 
simulation shows that the burning happens in the lam- 
inar regime at the beginning of flame acceleration, and 
considerable turbulence is generated close to the end of 
the process. Turbulence level depends typically on the 
obstacle spacing. 

We took slip and adiabatic boundary conditions at the 
tube wall: 



n • u = 0, n • VT = 0, 



(27) 



where n is the unit normal vector at the wall. At the 
open end of the channel, non-reflecting boundary condi- 
tions were used. As initial conditions, we used a hemi- 
sphere of hot "burnt" gas at the channel axis at the 
closed end of the tube, with the temperature profile 
given by the analytical solution of Zel'dovich and Frank- 
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FIG. 7: Scaled velocity of the flame tip versus the scaled 
tip position for planar and axisymmetri c g eo met ry. The the- 
oretical lines correspond to Eqs. |l6| , (24). The sym- 
bols show results of numerical modeling for a — 1/2; 2/3 
and Az/R = 1/4; 2. The experimental data shows results of 
Ref. [TQ. 



Kamenetskii [21125]. 



T f + (T b -T f )exp 



V 'x 2 + z 1 



L 



if z 2 



■ x < r 



f 



T 
Y 



2 ^ 2 

x >r /5 



67} if z l 
(T b -T)/(T b -T f ), P = P f , u x 



0, u z = 0. (28) 



Here rj is the radius of the hemisphere. The boundary 
of the hot gas is not a flame front yet, and we take both 
cold and hot gas initially at rest. The finite initial radius 
Tf is equivalent to a time shift, which requires proper 
adjustments when comparing the theory and numerical 
simulations. When necessary, we shifted the numerical 
solution in time to have the theory and the modeling 
results starting at the same point. 

A 2D hydrodynamic Navier-Stokes code adapted for 
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FIG. 8: Snapshots of temperature and velocity in the flow 
generated by an accelerating flame in an axisymmetric geom- 
etry for Ma = 0.001, Az/R = 1/4; a = 1/2. 



FIG. 9: Ratio of the acceleration rate obtained numerically 
to predicted theoretically, Eq. (7), versus the initial Mach 
number for the planar geometry and a = 1/3; 1/2; 2/3 and 
Az/R =1/4. 



parallel computation [29-31 was used. The numerical 
scheme is second-order accurate in time, fourth order in 
space for the convective terms, and second order in space 
for the diffusive terms. The code is robust and accu- 
rate, having been successfully used in combustion and 
aero- acoustic applications. The code is available in 2D 
(Cartesian and cylindrical axisymmetric) and 3D Carte- 
sian versions. In the present work, we performed only 2D 
simulations to save computational time and to be able to 
perform a large number of simulation runs required for a 
thorough investigation of the problem. 

A uniform grid with quadratic cells of size 0.2Lf was 
used to ensure isotropic propagation of the curved flame 
in x and y directions. The longitudinal size of the calcu- 
lation domain changes dynamically, following the leading 
pressure wave. Splines of the third order were used for 
re-interpolation of the flow variables during periodic grid 
reconstruction to preserve second-order accuracy of the 
numerical scheme. We performed several test simulation 
runs with resolutions of 0.125L/, 0.25L/, 0.5L/. The 
test demonstrated that the flame velocity grows expo- 
nentially for all chosen resolutions, with the difference in 
the scaled acceleration rate a not exceeding 6%. The de- 
pendence of the scaled flame tip velocity U t i p /Uf on the 
scaled distance Z tip /R is shown in Fig. [3] for different 
resolutions. Resolution tests also showed convergence of 
the acceleration rate a value with increasing resolution. 



4. SIMULATION RESULTS AND DISCUSSION 

Simulation of the flame acceleration in channels/tubes 
with obstacles was performed at various flow param- 
eters. In particular, we used various blockage ratios 
a = 1/3; 1/2; 2/3, various spacings between the ob- 
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FIG. 10: Scaled velocity of the flame tip versus tip position 
for the planar geometry and a = 1/3; 1/2; 2/3, Az/R =1/4, 
and different reaction order with respect to density n — 1,2. 
The plot with n = 2 shows also transition to detonation. 

stacles, Az/R =1/4; 1/2; 1; 2, and various initial Mach 
numbers. Figure [4] shows the characteristic flame shape 
and flow velocity obtained at various conditions. It is 
found that, in spite of these differences, all snapshots 
demonstrate the same basic feature of flame accelera- 
tion in channels with obstacles described theoretically in 
Sec. 2, namely, the flame tip propagates fast along the 
unobstructed part of the channel, leaving pockets of un- 
burned fuel mixture in between the obstacles. Delayed 
burning between the obstacles involves various levels of 
turbulence created by the flow. In Fig. [4] (a) the flow is 
laminar, and this regime is quite similar to the theoret- 
ical one sketched in Fig. [I] and described in Sec. 2. In 
contrast, Fig. [4] (g) shows rather strong turbulence cor- 
responding to the developed stage of flame acceleration. 
Figures [4] (c, e) show the flame shape for relatively large 
obstacle spacing: the snapshots are similar to the experi- 
mental photos of accelerating flames in Ref. [10 . We also 
observe a strong jet-flow developing in the unobstructed 
part of the channel in Fig. [4] (b, d, f, h) which is an im- 
portant part of the acceleration mechanism described in 
Sec. 2. Figure [5] shows the z- velocity component of the 
jet along the channel axis for the initial laminar stages of 
burning. As predicted by the theory, Eq. Q, plotted by 
the dashed line, we observe almost linear increase of the 
gas velocity from the closed end of the tube to the flame 
tip. The x-velocity component in Fig. [6] also demon- 
strates good agreement between the theory and simula- 
tion. Comparison of the theory and simulation in Figs. 
[5j [6] supports the potential flow model in the burnt gas 
employed in the analytical theory of Sec. 2 and of Ref. 
[T7] . The slight deviation between the theoretical and 
numerical results is due to the finite flame thickness and 
viscosity employed in the simulation. For example, be- 
cause of viscous friction, the z- velocity component in the 
simulation inevitably decreases at the border between the 




FIG. 11: Temperature field for three consecutive moments 
during deflagration-to-detonation transition for the second- 
order reaction for 8 = 8, Ma = 0.005, Az/R =1/4, a = 1/2. 



obstructed and unobstructed parts of the channel, which 
produces vorticity in that region as shown in Fig. 1 (c) of 
Ref. [17 . Another interesting question is related to the 
vorticity generated at the curved flame front as predicted 
by the classical theory [2] . In the present case the curved 
shape of the flame tip plays a minor role in the flame 
dynamics in comparison with the powerful acceleration 
mechanism described in Sec. 2. For example, replacing 
the curved flame tip by a planar one in Fig. [4] (a), we 
obtained negligible change of the total flame surface area 
as well as the total burning rate. Similarly, the vorticity 
generated by the curved flame tip plays a minor role in 
comparison with the strong jet flow described in Sec. 2 
and with the vorticity generated in that flow because 
of viscous forces. In this sense the present problem is 
completely different from, say, the Darrieus-Landau in- 
stability for which a curved flame shape and the flame 
velocity increase are intrinsically related to the vorticity 
production [32] . 

One of the most interesting questions in this study is 
to identify how the flame acceleration depends on the ob- 
stacle parameters: blockage ratio and spacing, i.e. pocket 
depth and width. According to the theory of Sec. 2, the 
acceleration rate is determined by the size of the unob- 
structed part of the channel, R(l — a). In the theoretical 
model, the pocket depth and width do not influence the 
acceleration, provided that the depth is not too small and 
the width is not too large, for which the notion of a pocket 
becomes meaningless. Taking proper scaling of the tip 
position Zf/R(l — a), we observe that the numerical re- 
sults reproduce the theoretical predictions for different 
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FIG. 12: Temperature profiles along the channel axis at 
different time instants close to the instant of detonation trig- 
gering; other parameters are the same as for Fig. |11| 



FIG. 13: Temperature profiles along the channel axis, at the 
obstacle edges and at the wall for Uft/R = 0.30661; other 
parameters are the same as for Fig. |11| 



obstacle parameters, see Fig. [7] Additional numerical 
data, which is not presented in Fig. [7] to avoid clutter, 
can be found in Ref. [17]. In the event of small spac- 
ing between the obstacles (Az/R = 1/4), the flame tip 
accelerates monotonically as described by the theoretical 
model. When the spacing is considerable (Az/R = 2), 
we observe strong turbulent pulsations with space pe- 
riod well-correlated with the distance between the ob- 
stacles. Still, even in this case, the average velocity of 
the flame tip is predicted by the theory of Sec. 2 with 
good accuracy. Figure [7] also shows good agreement of 
the theory and simulation with the experimental results 
[IT], which also demonstrate noticeable velocity pulsa- 
tions because of the large spacing between the obstacles. 
When compared to the obstacle positions (not shown in 
the figure because of the chosen scaling), both simulation 
and experiment demonstrate an increase in the flame tip 
velocity at every obstacle. Maximum pulsation velocity 
corresponds approximately to the middle between two 
obstacles. 

For the axisymmetric geometry, the theory of Sec. 2 
predicts considerably stronger flame acceleration as com- 
pared to the planar case. Figure [7] shows the veloc- 
ity of the flame tip versus the position in the tube for 
Ma = 10 -3 , a = 1/2; 2/3, and demonstrates much faster 
flame acceleration. The growth rate is about twice larger 
in agreement with Eqs. ([16]), (18). However, quantitative 



agreement between the theory and simulation in the ax- 
isymmetric case is not as good as in the planar geometry. 
In the planar case of Fig. [7] the quantitative difference is 
less than 5%, while in the axisymmetric case this differ- 
ence increases from about 10% at the beginning to 20% 



at later time. Equation (18) provides better agreement 



with the simulation than Eq. (25), which may indicate 



a minor role of the radius growth for the flame skirt in 
the pockets. This effect may be observed directly for 
the flame shape in the axisymmetric geometry shown in 



Fig. [8j which shows that even minor penetration of the 
flame skirt into the pockets between the obstacles pro- 
vides quite strong flame acceleration in the axisymmet- 
ric geometry. For comparison, the flame skirt penetrates 
much deeper into the pockets in the planar case of Fig. 
[4] (a). We also observe that the difference between theory 
and modeling becomes larger in Fig. [7] as the flame ve- 
locity increases. The quantitative difference of 10-20% in 
the axisymmetric case is partly related to the influence of 
viscosity and the moderate tube width, corresponding to 
moderate values of the Reynolds number. As obtained 
in numerical studies of the early acceleration of finger 
flames [23] , moderate tube width reduces flame accelera- 
tion as compared to the theory. However, there are also 
other reasons that cause a reduction of the flame acceler- 
ation in the numerical modeling. The first is the incom- 
pressibility assumption of the theory. This assumption 
is acceptable at the beginning of the acceleration, but 
deteriorates as the local Mach number increases as the 
flame accelerates. The deviation develops faster for the 
axisymmetric geometry because of its faster acceleration. 
In order to study the influence of gas compression on the 
flame acceleration, we investigated the dependence of the 
scaled acceleration rate a on the initial Mach number, as 
shown in Fig. [9] for various values of the blockage ratio. 
It is thus seen that the scaled acceleration rate decreases 
strongly with the Mach number. For initial Mach num- 
ber Ma = 10 -2 instead of Ma = 10 -3 , the acceleration 
rate is approximately reduced by a factor of two. 

The slowdown of flame acceleration because of gas 
compression agrees with the concept that the flame prop- 
agation velocity cannot exceed the value of the Chapman- 
Jouguet (CJ) limiting state, for which the downstream 
flow is sonic. We therefore expect saturation of the flame 
tip velocity to a certain steady value at the end of the ac- 
celeration process, but prior to an explosion. Indeed, Fig. 



10 demonstrates for the planar geometry such a satura- 
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tion, obtained in the numerical simulation at the final 
stage of flame acceleration, for various blockage ratios. 
The initial Mach number in Fig. [To] is Ma = 5 • 1CT 3 , 
and the saturation velocity is about (2.5 — 3.0)c s . For 
comparison, one-dimensional theory predicts the CJ de- 
flagration velocity with respect to the wall in the limit of 
large energy release (6 >> 1) as [T8l [T9] 



U C j 




(29) 



In the present case of 7 = 1.4, = 8, Eq. (29) yields 
Ucj ~ 2.68c s , which is about the saturation velocity ob- 
tained in the present numerical simulation. More accu- 
rate theoretical calculation of the C J deflagration velocity 
without the assumption of O >> 1 yields Ucj ~ 3.38c s , 
which is also close to the present numerical results. 

Finally, we discuss how flame acceleration in channels 
with obstacles may lead to DDT. Since the time and po- 
sition of DDT are quite sensitive to the chemical kinetics 
adopted, our results on detonation triggering should be 
considered as qualitative. At the same time, they demon- 
strate the general features of DDT irrespective of the 
particular fuel mixture. It is well known that any flame 
propagating from a closed end pushes a flow in the fuel 
mixture with a weak shock/compression wave at the head 
of the flow. The flame acceleration renders the compres- 
sion wave stronger, until it develops into a shock of con- 
siderable amplitude. Preheating of the fuel mixture by 
the shock is conventionally considered as one of the main 
elements of DDT both in smooth tubes and in channels 
with obstacles [TUB]. The temperature behind the shock 
increases and the reaction time in any compressed gas 
parcel decreases drastically. The decrease in the reaction 
time may result in explosion and DDT ahead of the flame 
front unless the parcel is burnt by the flame before active 
explosion is initiated. Thus, in general, we may expect 
two possible outcomes for the flame acceleration: 1) If 
the reaction time behind the shock is sufficiently short, 
then it drives the explosion and DDT; 2) The reaction 
time may be longer than the interval available for a gas 
parcel to travel between the shock and the flame. In this 
case explosion does not occur and the final state of flame 
acceleration is the CJ deflagration. For comparison, the 
possibility of spontaneous explosion ahead of an acceler- 
ating flame was considered in the theory [33] for smooth 
tubes. Both CJ detonation and deflagration have also 
been found in smooth tubes experimentally in Ref. [16]. 
In channels with obstacles, the state of CJ deflagration 
is also known as " fast flames" [6j [8] . In the present sim- 
ulations we observed both possibilities of DDT and CJ 
deflagration for different reaction kinetics. Taking reac- 
tion of the first order with respect to density (designated 
by n = 1 in Fig. 10), we obtained statistically steady 
CJ deflagration at the end of flame acceleration with no 
explosion or DDT, see Fig. 
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10 This result indicates that 



(c) 



FIG. 14: (a) Temperature, (b) pressure gradient modulus and 
(c) pressure fields close to the point of explosion triggering for 
the second-order reaction for B = 8, Ma = 0.005, Az/R = 
1/4, a = 1/2. 



the decrease in the reaction time behind the shock is 
not sufficient, and gas parcels are consumed by the flame 
front before spontaneous reaction develops into a power- 
ful explosion. Thus, in order to observe DDT, we need 
to take another reaction mechanism, which is more sen- 
sitive to pressure and temperature increase in the shock. 
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Similar to Ref. [26], we considered a reaction of the sec- 
ond order with respect to density, n = 2, and obtained 
explosion triggering and DDT, see Fig. [Toj Remarkably, 
in this case the reaction rate is so sensitive to pressure 
and temperature that the DDT occurs before the flame 



reaches the CJ deflagration state. Figure 11 shows char- 
acteristic temperature snapshots at the DDT. It is seen 
that the accelerating flame acts like a piston pushing a 
shock, which raises the temperature of the fuel mixture 

(a, b). The 
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ahead of the flame, as shown in Fig. 
snapshot of Fig. [IT] (c) already corresponds to detona- 
tion. Development of the temperature profiles ahead of 
the flame front along the channel axis is shown in Fig. 



12] for the time instants Uft/R = 0.30613 - 0.30844 close 
to the DDT time. Unlike the case of smooth tube with 
monotonic temperature distribution in the compression 
wave [33 , in Fig. 12 we observe noticeable tempera- 



ture pulsation due to secondary shocks reflected from the 
obstacles. The main tendency nevertheless remains the 
same: we can see a considerable temperature jump in the 
main shock followed by further temperature increase in 
the compression wave from the shock to the flame front. 
The compression wave and the shock become stronger as 
the flame accelerates, until explosion starts and develops 
into detonation. 

Recent papers [20j [34] on DDT in tubes with smooth 
adiabatic wall have also demonstrated the important role 
of viscous heating at the wall in addition to shock heat- 
ing. Because of viscous heating, the temperature at the 
wall of a smooth adiabatic channel is larger than that at 
the axis, and numerical modeling demonstrates DDT on- 
sets at the wall. A similar physical mechanism of viscous 
heating may also be identified in channels with obstacles, 
though with proper modifications due to the specific ge- 
ometry. To elucidate the mechanism, Fig. [13] shows the 
temperature profiles for Uft/R = 0.30661 along the chan- 
nel axis, along the wall and at specific obstacle edges. 
We can see that on average, temperature is considerably 
higher at the channel axis. Obstacles reduce the shock 
strength, which results in much lower average temper- 
ature at the obstacle edges, with the lowest tempera- 
ture at the wall deep in the pockets. Still, obstacles not 
only moderate the main shock, but they also produce 
hot spots, which may be crucial in explosion triggering 
and DDT [2j [35]. Though the average temperature is 
lower at the obstacle edges, we also observe sharp peaks 
of temperature ahead of every obstacle, which are much 
higher than the respective temperature at the channel 
axis. There are two possible reasons for these tempera- 
ture peaks. First, secondary shock waves reflected from 
the obstacles may produce local temperature increase. 
Second, the strong jet-flow pushed by an accelerating 
flame slows down at the obstacles and generates vortices 
in the pockets with high velocity gradients. Slowdown 
of the jet flow increases local pressure and temperature. 
Viscous dissipation of the vortices leads to additional 



temperature increase, which has the same effect as vis- 
cous heating at the wall in smooth tubes [20] |34] . The 
important role of viscous heating in producing tempera- 
ture peaks is especially obvious in Fig. [13] at the positions 
Z/R = 25.75; 26; 26.25, which are ahead of the strong 
shock position, Z/R = 25.44. At this simulation run ex- 
plosion starts at an obstacle edge, as shown in Fig. [14 
It is noted the exact position of explosion triggering may 
depend on the particular obstacle geometry, e.g. on the 



spacing between the obstacles. In Figs. 11 , 14 we used a 
rather small spacing Az/R =1/4. Other simulations [12] 
performed for larger spacings also demonstrated the pos- 
sibility of explosion triggering deep in the pockets with 
the hot spots produced by reflected shocks. Thus, both 
in Ref. [12] and in the present simulation, obstacles play 
an important role in explosion triggering and DDT. 



5. SUMMARY 

This paper presents the theory and numerical simu- 
lation of flame acceleration in channels with obstacles. 
We showed theoretically as well as computationally that 
flame acceleration is noticeably stronger in the axisym- 
metric geometry as compared to the planar one. We also 
considered the influence of gas compression on the flame 
acceleration, and showed numerically that the flame ac- 
celeration rate decreases with increasing initial Mach 
number, and that the velocity of the accelerating flame 
eventually saturates to a value that is supersonic with 
respect to the wall and is correlated to the known CJ 
deflagration speed. This saturation state has been re- 
ferred to as that of fast flames in experimental studies 
[6j [21] . We also demonstrated numerically the possibility 
of DDT in the geometry of obstructed channels. 
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